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I .  INTRODUCTION 

In  some  recent  work,*  we  exposed  a  one-dimensional  model  and  developed 
a  formalism  for  studying  plasma  dynamics  in  an  arc-driven  rail  gun.  A  set  of 
governing  equations  was  derived  and  solved  numerically  to  determine  the  prop¬ 
erties  of  the  arc  for  various  rail  guns,  projectile  masses,  and  accelerating 
characteristics.  Approximate  analytic  solutions  to  the  resulting  equations 
were  also  obtained  and,  from  these  solutions,  some  scaling  relations  were  de¬ 
rived.  The  extent  to  which  the  resulting  solutions  accurately  represent  the 
arc's  characteristics  was  assessed,  and  the  necessity  for  extending  the  work 
to  two  dimensions  was  emphasized.  In  this  report  that  extension  is  under¬ 
taken. 

The  rationale  for  the  Army's  interest  in  rail  guns  was  thoroughly  dis¬ 
cussed  in  Ref.  1  and  will  not  be  discussed  here.  It  is  well  to  point  out, 
however,  that  the  most  successful  rail-gun  experiments  undertaken  thus  far 
have  employed  a  plasma  arc  to  drive  the  projectile.  No  doubt  the  success  of 
these  experiments  results  from  the  ability  of  the  arc  to  maintain  better  con¬ 
tact  with  the  rails  than  can  the  sliding  contacts  used  in  solid-armature  rail 
guns.  At  the  higher  velocities  anticipated  in  future  experiments  an  arc  arma¬ 
ture  is  likely  to  be  indispensable.  Consequently,  it  is  important  to  assess 
the  properties  of  the  arc,  to  determine  its  erosive  effects  on  the  rails  and 
projectile,  and  indeed  to  determine  the  feasibility  of  using  arcs  in  rail  guns 
sufficiently  large  to  be  of  military  interest. 

Before  turning  to  the  calculations  in  this  report  we  will  summarize  some 
of  the  more  recent  developments  in  rail-gun  technology.  Westinghouse,  under 
DARPA/ARRADCOM  support,  has  nearly  completed  construction  of  a  large  gun  in¬ 
tended  to  accelerate  a  300-g  projectile  to  3  km/s  in  about  3  m.  The  current 
intention  in  that  experiment  is  to  use  a  solid  armature  in  the  initial  stages 
but,  after  transfer  of  the  gun  to  Dover,  some  experiments  employing  an  arc  are 
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intended.  Arc-driven  rail-gun  experiments  have  been  undertaken  most  recently 
by  a  team  from  Los  Alamos  and  Livermore.  They  have  constructed  rail  guns  which 
use  a  magnetic-flux  compression  generator  as  a  primary  power  supply,  and  have 
successfully  accelerated  projectiles  of  a  few  grams  to  several  kilometers  per 
second.  In  addition,  in  one  experiment,  a  165-g  projectile  was  accelerated  to 
about  350  m/s.  The  latter  effort  represents  the  most  massive  projectile  yet 
accelerated  by  a  rail  gun.  In  addition  to  these  experiments,  BRL  has  recently 
initiated  an  experimental  program  and  has  constructed  a  small  rail  gun  powered 
by  a  capacitor  bank.  The  intention  of  the  present  and  near-future  experiments 
at  BRL  will  be  to  analyze  the  properties  of  the  arc  in  some  detail,  rather 
than  to  accelerate  ever-more-massive  projectiles  to  ever-higher  velocities. 

1.  John~D'.  Powell  and  Jad  H.  Batteh ,  ", Plasma  Dynamics  of  the  Arc- Driven  Rail 
Gun , "  Ballistic  Research  Laboratory  Report  No.  ARBRL-TR-02267 ,  September 
1980.  See  also  "Plasma  Dynamioe  of  an  Arc-Driven,  Electromagnetic,  Pro¬ 
jectile  Accelerator,"  J.  Appl.  Phye.  52,  2717  (1981).  (AP  A092345) 

2.  R.S.  Hawke,  A.L.  Brooks,  F.J.  Deadriak,  J.K.  Scudder ,  C.M.  Fowler,  R.S. 
Caird ,  and  D.R.  Peterson,  "Results  of  Rail  Gun  Experiments  Powered  by 
Magnetic  Flux  Compression  Generators ,"  IEEE  Trans.  Magnetics  18j  821 
(1982). 
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Some  diagnostic  measurements  have  been  undertaken  already,  and  more  sophisti¬ 
cated  ones  are  planned  for  the  immediate  future. 

In  the  present  calculation  we  wish  to  employ  the  same  model  and  assump¬ 
tions  that  were  used  in  Ref.  1.  A  schematic  diagram  of  the  model  is  shown  in 
Fig.  1.  Sides  one  and  two  represent  the  rails,  which  are  infinitely  extended 
in  the  z  direction,  and  which  carry  current  in  the  direction  indicated.  Side 
three  represents  the  power  supply  and  side  four  contains  both  the  arc  (between 
x  =  xQ  and  x  =  i)  and  the  projectile  (shaded  area).  The  resulting  current 

configuration  produces  a  magnetic  field  in  the  space  between  the  rails,  and 
this  field  interacts  with  the  current  through  the  arc,  accelerating  both  it 
and  the  projectile  in  the  x  direction. 

The  intention  of  the  calculation  is  to  determine  the  fluid-dynamical  pro¬ 
perties  of  the  arc,  the  electromagnetic  fields  associated  with  it,  and  the 
acceleration  imparted  to  the  arc  and  projectile.  In  Ref.  1,  these  properties 
were  assumed  to  vary  only  in  the  direction  of  acceleration,  while  in  the  pre¬ 
sent  work  variations  in  one  transverse  direction,  the  y  direction,  will  also 
be  allowed.  The  additional  assumptions  employed  and  justified  in  Ref.  1  will 
also  apply  to  this  calculation:  First,  since  the  conductivity  of  the  rails  is 
several  orders  of  magnitude  higher  than  that  of  the  arc,  we  simplify  the  prob¬ 
lem  by  assuming  that  the  rails  are  actually  perfect  conductors.  The  finite 
conductivity  of  the  arc,  however,  is  accounted  for.  Second,  we  assume  that 
the  plasma  is  at  most  doubly  ionized  and  neglect  higher  degrees  of  ionization. 
This  assumption,  while  not  fundamental  to  the  calculation,  greatly  simplifies 
the  algebra  and  is  probably  not  unreasonable  in  the  temperature  ranges  under 
consideration.  Third,  we  assume  that  energy  losses  within  the  arc  occur  only 
by  radiation  and  neglect  ordinary  heat  conduction.  As  was  demonstrated  ex¬ 
plicitly  in  Ref.  1,  heat  conduction  is  negligible  at  the  temperatures  which 
characterize  arcs  in  rail  guns  of  experimental  interest.  Fourth,  we  assume 
that  the  current  per  unit  height  on  the  surface  of  the  rails  is  constant  in 
time.  This  condition  can  be  satisfied  in  practice  by  using  as  a  power  source 
an  inductive  store  having  a  sufficiently  high  inductance.  Wo  will  also  assume 
in  the  present  treatment  that  the  gas  is  inviscid. 

In  addition  to  the  above  assumptions,  we  also  assumed  in  the  one¬ 
dimensional  calculations  that  the  properties  of  the  arc  were  steady  in  a  frame 
of  reference  that  accelerated  with  the  projectile,  and  demonstrated  that  the 
governing  equations  possessed  a  solution  in  that  limit.  In  the  two-dimensional 
case,  as  we  shall  see,  no  such  solution  exists  for  the  steady  situation  in 
which  the  arc  is  at  rest  relative  to  the  projectile.  However,  if  the  effect 
of  the  arc's  acceleration  upon  its  fluid-dynamical  properties  can  be  neglected, 
then  a  steady  solution  similar  to  the  one-dimensional  case  does  exist.  We 
shall  assume  this  to  be  the  case  in  our  actual  numerical  calculations,  al¬ 
though  we  will  derive  a  time- dependent  formalism  in  order  to  make  possible 
more  general  calculations  in  the  future.  The  physical  reasons  for  the  differ¬ 
ences  between  the  one-  and  two-dimensional  cases  are  discussed  in  Sec.  Ill  of 
the  report,  and  some  discussion  of  the  significance  of  neglecting  the  accel¬ 
eration  is  given  in  Sec.  IV. 

Before  discussing  the  present  work,  it  is  necessary  to  indicate  the  ra¬ 
tionale  for  undertaking  a  two-dimensional  calculation  and  what  differences 
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might  be  expected  with  regard  to  the  one- dimensional  results.  In  Ref.  1  it 
was  demonstrated  that  of  all  the  energy  deposited  into  the  model  rail-gun  sys¬ 
tem  by  the  power  supply,  slightly  less  than  half  of  the  onergy  became  asso¬ 
ciated  with  kinetic  energy  of  the  arc  and  projectile.  An  equal  amount  was 
stored  in  the  magnetic  field  and  the  remaining  onergy,  a  small  fraction  of 
the  total,  was  dissipated  in  the  arc.  In  a  two-dimensional  calculation  a 
similar  conclusion  should  hold  and,  since  energy  must  still  be  conserved,  the 
projectile  will  have  essentially  the  same  kinetic  energy  and  an  equivalent 
amount  of  energy  should  be  stored  in  the  field.  Consequently,  the  arc  pres¬ 
sure,  which  accelerates  the  projectile,  will  have  to  be  similar.  On  the  other 
hand,  it  is  evident  that  more  surface  area  is  available  for  radiation  in  the 
two-dimensional  calculation,  so  the  temperature  of  the  arc  should  be  signifi¬ 
cantly  smaller.  Similarly,  the  mass  density  will  be  larger  in  the  two- 
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Figure  1.  Model  for  Arc  Calculations 
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dimensional  case,  since  it  must  increase  as  the  temperature  decreases  to  satis¬ 
fy  the  equation  of  state  (assuming  constant  pressure  and  degree  of  ionization). 
It  is  clear,  therefore,  that  a  one- dimensional  model  probably  overestimates 
the  temperature  in  the  arc,  underestimates  the  density,  and  gives  fairly  ac¬ 
curate  results  for  the  pressure  and  fields.  However,  a  two-dimensional  model 
not  only  allows  for  radiation  at  the  front  and  the  back  of  the  arc,  but  also 
along  the  surfaces  in  contact  with  the  rails.  Since  these  latter  surfaces  are 
usually  large  in  experiments  of  interest,  a  significant  improvement  in  the 
temperature  and  density  calculations  should  be  anticipated. 


In  addition,  it  was  found  that  the  one-dimensional  calculation  predicted 
a  rather  anomalous,  nonphysical  manner  in  which  the  temperature  varied  with 
the  length  of  the  arc.  Intuitively,  one  expects  that  as  the  length  increases 
the  temperature  should  fall,  since  the  resistance  of  the  arc  becomes  smaller. 
However,  it  was  shown  in  Ref.  1  that  the  temperature  remained  essentially  con¬ 
stant  as  the  arc  length  was  increased.  The  reason  for  this  unexpected  behavior 
was  that,  as  the  length  increased,  the  temperature  gradients  within  the  arc 
also  became  smaller.  Consequently,  the  energy  dissipated  was  radiated  away  less 
efficiently.  In  the  one-dimensional  problem  the  two  effects  (smaller  dissipa¬ 
tion  and  less  efficient  radiation)  almost  exactly  compensated  each  other  with 
the  result  that  little  variation  in  the  temperature  was  observed.  As  was 
pointed  out,  however,  this  unusual  effect  was  just  a  peculiarity  of  the  one¬ 
dimensional  model  and,  for  a  two-  or  three-dimensional  case,  the  temperature 
should  be  expected  to  fall  with  increasing  arc  size.  The  present  calculation, 
therefore,  should  predict  in  some  detail  how  that  variation  occurs. 


The  organization  of  the  report  is  the  following.  In  Sec.  II  the  time- 
dependent  formalism  for  treating  the  two-dimensional  problem  is  developed  and 
the  appropriate  set  of  governing  equations  derived.  In  Sec,  III,  the  approxi¬ 
mations  employed  in  the  solution  of  the  equations  are  discussed,  and  numerical 
as  well  as  approximate  analytic  solutions  to  the  equations  are  obtained  for 
several  cases  of  interest.  In  addition,  a  set  of  approximate  scaling  relations 
is  derived  which  indicate  how  arc  properties  vary  with  gun  size,  projectile 
mass,  and  acceleration  characteristics.  Finally,  in  Sec.  IV,  some  discussion 
of  the  results  is  given. 


II.  MATHEMATICAL  FORMALISM 

In  this  section  we  will  be  concerned  with  deriving  the  equations  which  de¬ 
termine  the  electromagnetic  fields  associated  with  the  arc,  its  fluid-dynamical 
characteristics,  and  the  acceleration  of  the  arc  and  projectile.  A  general, 
time-dependent  formalism  will  be  worked  out  here,  though  additional  approxima¬ 
tions  will  be  made  in  the  following  section  where  the  solution  is  carried  out. 

Applying  Maxwell's  equations  in  the  space  between  x  =  x.  and  x  «  1  in 
Fig.  1,  we  have 


V 

X  B  «  yj 

(2.1) 

V 

c  3B 

x  E  -  -  at 

(2.2) 
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and 


J  *  o(E  +  v  x  B)  .  (2.3) 

In  these  equations  1*  represents  the  magnetic- induction  field,  $  the  electric- 
field  intensity,  J  the  current  density,  v  the  plasma  material  velocity,  o  the 
plasma  conductivity,  and  y  the  permeability  which  is  assumed  to  be  that  of 
free  space.  The  displacement  current  has  been  neglected  in  Eqs.  (2.1)-(2.3), 
since  the  propagation  velocities  under  consideration  are  small  compared  to 
the  speed  of  light. 

It  is  convenient  to  make  a  transformation  to  a  dimensionless  coordinate 
system  which  accelerates  with  the  arc.  Specifically,  we  let 

x-xQ 

5  -  — —  n  =  y/w  ,  (2.4) 

a 

where  Z  «  S.  -  x_  is  the  distance  from  the  projectile  to  the  trailing  edge  of 

fl  vJ 

the  arc.  The  transformation  is  then  effected  by  noting  that 


3__  =  1 _ a 

3x  ~  Z  35  3y  ""  w  3n 

a  (2.5) 

_3  J3_  ^0  __3  £_  ^a  J3 

3t  “  3t"  l  '  35  "  %  3t'  35 

a  a 

where  v^  is  the  velocity  of  the  trailing  edge  of  the  arc  and  where  the  nota¬ 
tion  3/3t'  is  used  to  denote  the  time  derivative  measured  in  the  moving  frame. 

We  now  take  the  curl  of  Eq.  (2.1),  note  that  $  •  ^  is  zero,  and  use  Eq. 
(2.2)  to  obtain  a  differential  equation  for  1$  in  terms  of  the  conductivity 
and  plasma  velocity.  After  considerable  simplification  we  find  that  the  ap¬ 
propriate  equation  in  the  new  coordinate  system  is 


.2  2  3B 

*a  “  at'  =  “ 


-  y  a  SL 

t 


w  w 


o  D  nL 

— J  +  *a 
35 
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(2.6) 


In  writing  Eq.  (2.6)  we  have  noted  that,  from  symmetry,  $  has  only  a  single 
component  in  the  z  direction  and  have  defined  v'  as  the  plasma  velocity  rela 
five  to  the  trailing  edge,  i.e. , 
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Once  the  induction  field  has  been  determined  from  Eq.  (2.6),  the  current  den¬ 
sity  J  and  electric-field  intensity  E  then  follow  directly  from  Eqs.  (2.1) 
and  (2.3),  respectively. 

The  fluid-dynamic  properties  of  the  arc  are  determined  largely  from  equa¬ 
tions  which  express  the  conservation  of  mass,  momentum,  and  energy  for  the 
arc  as  a  whole.  Specifically,  we  have 


3p 

at 


+ 


1  • 


(P  V)  .  0 


(2.8) 


P 


pv 


$  V  ♦  $  •  P  =  * 

'V 


(2.9) 


p  +  .  $  e  ♦  3  •  q  +  P:^v 

% 

In  these  equations  p  represents  the 
energy;  JP  is  the  second-order  pressure  or  momentum-flux  tensor;  q  is  the  heat 
flux  which  accounts  for  losses  due  to  radiation  and  conduction;  finally,  ?  is 
the  force  per  unit  volume  acting  on  the  plasma  and,  for  the  problem  at  hand, 
it  is  given  by 


=  J 2  fa 


(2.10) 


plasma  density  and  e  its  specific  internal 


*  =  J  xl 


(2.11) 


As  before,  it  is  convenient  to  transform  Eqs.  (2.8)-(2  10)  to  the  moving 
frame.  Using  Eq.  (2.5),  we  find 
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(2.13) 


Unprimed  quantities  will  be  used  throughout  to  refer  to  quantities  measured 
in  the  fixed  frame  and  primed  variables  to  those  measured  in  the  moving  frame. 
In  the  nonrelativistio  limit  we  have,  for  example, 
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and 
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(2.15) 


In  obtaining  the  above  results  we  have  replaced  £  by  the  hydrostatic  pressure 

P  since  the  gas  in  inviscid,  have  used  Eqs.  (2.1)  and  (2.11)  to  express  the 
force  F  solely  in  terms  of  the  induction  field,  and  have  denoted  by  a  the 
acceleration  of  the  trailing  edge  of  the  arc.  The  resulting  equations  are 
identical  to  those  generally  encountered  in  fluid  mechanics  except  for  the 
forces  which  are  electromagnetic  in  origin  and  except  the  term  proportional 

to  J  /o  which  accounts  for  the  Joule  heating  of  the  plasma.  It  should  also 
be  noted  that  the  term  which  varies  with  a  is  an  effective  force  on  the  plasma 
that  arises  because  the  reference  frame  is  not  inertial;  such  a  force  does 
not  arise  in  the  more  customary  laboratory- frame  calculations. 


Equations  (2. 12)- (2. 15)  are  insufficient  to  determine  the  flow  parameters 
uniquely  and  we  need,  in  addition,  the  internal  energy  as  a  function  of  P  and 
the  temperature  T,  an  equation  o£  state  relating  P  to  T  and  p,  and  finally  a 
representation  of  the  heat  flux  q.  For  simplicity,  we  asstime  an  ideal-gas 
equation  of  state  so  that 


P 


(1  *  Xx  +  2  x2>  kB  T  p 
m0 


(2.16) 


where  n>0  is  the  atomic  mass  of  the  ions  or  neutrals  constituting  the  arc,  kR 

is  Boltzmann's  constant,  and  x^  and  x^  are  th®  concentrations  of  singly  and 

doubly  ionized  atoms,  respectively.  Thtsse  concentrations  can  be  determined 
from  the  Saha  equations  which,  for  a  singly  and  doubly  ionized  gas,  can  be 
3 

written 


3.  Y.B.  ZeV dovish  and  Y.P.  Paizer ,  Physios  of  Shook  Waves  and  High-Temperature 
Hydrodynamic  Phenomena ,  ( Academic ,  New  York,  1966),  Vol.  I,  Chap.  3. 
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x.  (1  ♦  x.  +  2x^ 

1  1  4* 


*  K2  (T,P) 


(2.17) 


where 


,  Z.  ,  m  3/2  5/2  -I.  .  _ 

VF-z^r1^  (kBn'  8  1  *  ‘ 

1J 


(2.18) 


In  Eq.  (2.18),  m  is  the  electron  mass,  h  is  Planck’s  constant  divided  by  2ir, 
and  1^  is  the  energy  needed  to  ionize  the  atom  j  times.  The  quantity  re¬ 
presents  the  electronic  partition  function  for  the  jth  ior.;  specifically,  it 
is  given  by 


Z.  =  J  g..  e"Uji/kBT 

J  v  J1 


(2.19) 


where  u^  is  the  energy  of  the  ith  electronic  state  of  the  jitt  ion  and  where 

g..  is  the  appropriate  degeneracy  factor  for  this  level .  Values  of  these 
J1  4  1 

parameters  for  copper  and  its  first  two  ions  have  been  tabulated  elsewhere.  * 

In  the  temperature  range  with  which  we  will  be  concerned,  the  major  source 
of  energy  loss  within  the  arc  is  radiation  and  ordinary  heat  conduction  can  be 
neglected.  Furthermore,  if  the  radiation  mean  free  path  is  small  compared  to 
the  dimensions  of  the  arc,  as  we  will  find  to  be  the  case,  the  radiation  flux 

5 

can  be  represented  by 


q  =  - 


16  og  XT' 


(2.20) 


where  oc  is  Stefan’ s  constant  and  X  is  the  mean  free  path.  Expressions  for  X 
5  ^ 

have  been  worked  out  and,  for  the  problem  at  hand,  we  will  use  the  relation 


X(C,n)  = 


0.91  x  1011  kgT3  (l+Xj+Zxp 
P[l-x1-x2Te"Il^BT  4  4x1e_I2/kBT  +  9x2  e"  V*l?] 


. (2.21) 


T.  C.  E.  Moore t  Atomic?  Energy  Levels t"  National  Bureau  of  Standards  Ciro.  No. 
467 t  Vol.  II j  Washington,  DC,  1952. 

5.  See  Ref.  3,  Chap.  2. 

6.  Y. P.  Raizer,  "Simple  Method  for  Computing  the  Mean  Range  of  Radiation  in 
Ionized  Gaeee  at  High  Temperatures,"  Sov.  Phys.-JETP  37j  769  (I960).  See 
also  Ref.  3,  Chap.  5. 
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The  remaining  information  needed  to  complete  the  system  of  equations  des¬ 
cribing  the  arc  properties  8Te  expressions  for  the  arc  length  l  and  the  con- 

& 

ductivity  o.  If  we  define  as  the  mean  plasma  density  multiplied  by  its 


length  along  the  x  axis,  we  have  the  obvious  condition 

o. 


*a  “ 


a 


/ d$  / dn  p  (£ ,  n) 


(2.22) 


For  all  cases  of  interest  in  the  present  problem,  we  will  find  that  the  plasma 
is  almost  completely  ionized  with  few  neutrals  being  present  in  the  arc.  Con¬ 
sequently,  the  conductivity  can  be  well  approximated  by  that  for  a  completely 


7-9 

ionized  gas,  namely, 

2.63  x  10"2yp 
o(S,n)  *  - 2 - 

where  Z,  given  by 


/l , 23  x  107  T3/2\ 

v  / 


-i 


(2.23) 


Z  « 


X1  +  4x2 
________ 


(2.24) 


is  a  position-dependent  number  lying  between  one  and  two  which  is  indicative 
of  the  degree  of  ionization  in  the  gas  and  can  be  approximated  by  the  value 

0.6833.  Finally  the  parameter  n  is  the  electron  density  in  the  gas,  viz., 

e 


n  = 
e 


p(Xj  + 


2x2) 


m. 


(2.25) 


Equations  (2.1),  (2.3),  (2.6),  (2.12),  (2.13),  (2.14),  (2.15),  (2.16),  (2.17) 
and  (2.22),  together  with  t.he  supplementary  information  contained  in  this  sec¬ 
tion,  should  be  sufficient  to  determine  the  unknowns  J,  E,  B,  p,  v  ',  v  T, 

x  y 

p,  Xj,  x^  and  t,  .  In  the  following  section  several  solutions  to  these  equations 
are  carried  out. 


7.  L.  Spitzer,  Pkysiae  of  Fully  Ionized  Gaea 8,  ( Inters aience,  New  York ,  1965), 
Chap.  5. 

8.  R.S.  Cohen,  L.  Spitzer ,  and  P.  McR.  Routly ,  "The  Electrical  Conductivity  of 
an  Ionized  Gas,"  Rhys.  Rev.  80j  230  (1950). 

9.  L.  Spitzer  and  R.  Harm,  "Transport  Phenomena  in  a  Completely  Ionized  Gae , " 
Phye.  Rev.  89, ,  977  (1953). 
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III.  SOLUTION  OF  EQUATIONS 
A.  General  Remarks  and  Additional  Approximations 


In  Ref.  1  we  showed  that  for  a  one-dimensional  case  (no  dependence  on  y) 
the  equations  of  the  preceding  section  possessed  a  steady-state  solution  in 
which  v*  vanished,  and  we  justified  the  use  of  that  approximation  in  analysing 
the  arc's  properties.  If  we  make  a  similar  assumption  in  the  two-dimensional 
analysis,  Eqs.  (2.6)  and  (2. 13)- (2. 15)  become 


2  32B  „2  32B  2  3  log  a  33  ,2  3  log  a  3B 

w  +  *a  "  w  -“Tr-  W  +  *a  “Tn 


3P 

3P 

3n 


B  3B 

•  u  sr 

B  35 

y  W 


-  p  l&  a 


(3.1) 

(3.2) 

(3.3) 


and 


6qx  3qy  *> 

^+la?rslaw  J"/o 


(3.4) 


In  these  equations,  a  represents  the  common  acceleration  of  the  arc  and  pro¬ 
jectile,  namely, 


a  ~ 


u  j2 

2  to*.  *  P»  ’ 


(3.5) 


where  is  the  projectile  density  multiplied  by  its  length  along  the  x  axis. 
P 

Equation  (3.1)  is  to  be  solved  subject  to  the  boundary  conditions 


B(0,n)  =  yj 


Bci>n)  ■  °  (4-o  '(4,.i  ■  ° 


(3.6) 


The  conditions  on  B  are  easily  derived  from  Ampere's  law,  whereas  those  on 


3B 


are  necessary  in  order  that  the  x  component  of  the  current  density  vanish  at 
the  rail-arc  interface;  otherwise,  the  associated  electric  field  would  produce 
an  infinite  current  in  the  perfectly  conducting  rails.  The  conditions  obeyed 
by  P  and  T  are 


1C 


P(0, n)  -  0 

TCO.n)  -  TCI, n)  =  T(€,0)  *  T(e,l)  -  Tb  .  (3.7) 


The  latter  condition,  that  the  arc  is  in  contact  with  a  thermal  reservoir  at 
temperature  T^,  is  somewhat  different  from  the  boundary  condition  employed  in 

the  calculations  in  Ref.  1.  As  has  been  discussed  previously,  however,  the  re¬ 
sults  should  be  relatively  independent  of  boundary  temperature  and  the  latter 
condition  is  somewhat  easier  to  treat  in  the  analysis  which  follows.  This 
point  is  discussed  further  in  Sec.  Ill  C. 


If  we  now  integrate  Eq.  (3.2)  and  employ  Eqs.  (3.6)  and  (3.7),  we  find 


a 


r»  2 
/  o  d£  +  "• 

0 


(3.8) 


Substituting  Eq.  (3.8)  into  Eq.  (3.3)  leads  us  to  the  conclusion  that  a  steady 
solution  exists  only  if  p  is  independent  of  n  (one-dimensional  limit)  or  for 
vanishingly  small  values  of  the  acceleration  a.  The  physical  reason  that  no 
steady  solution  (with  the  gas  at  rest)  can  be  found  is  that  there  exists  an 
effective  force  on  the  arc  arising  from  the  last  term  on  the  right-hand  side 
of  Eq.  (3.2).  Although  this  force  lies  along  the  x  direction,  it  will  produce 
a  torque  on  the  arc  if  the  density  varies  in  the  y  direction  and  will  set  the 
arc  in  rotational  motion. 


Including  this  rotational  motion  of  the  arc,  as  well  as  possible  time- 
dependent  effects,  greatly  complicates  the  problem.  Therefore,  as  a  first  ap 
proximation,  we  will  neglect  the  effect  of  the  arc's  acceleration  upon  its 
fluid-dynamical  properties,  set  a  =  0  in  Eq.  (3.8),  and  find  the  steady  solu¬ 
tion.  Consequently,  the  solution  which  follows  corresponds  physically  to  a 
situation  in  which  the  projectile  is  held  fixed,  but  the  accelerating  current 
is  allowed  to  flow.  The  nature  of  the  approximation  is  discussed  in  greater 
detail  in  Sec.  IV. 

Making  use  of  the  approximation  discussed  above  then,  we  find  from  Eq. 

(3.8) 


and  this  result  satisfies  Eq.  (3.3)  automatically.  In  the  most  general  case, 
no  additional  simplification  of  the  equations  results.  However,  if  the  con¬ 
ductivity  can  be  approximated  by  a  separable  function  of  £  and  n  so  that 


c*(£,n)  =  a  (£)  a  (n)  , 
A  X 


(3.10) 


Eq.  (3.1)  can  be  solved  by  separation  of  variables. 


We  find 


B  «  y,j 
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°v  (C)  dC 


(3.11) 


Remarkably,  the  result  is  independent  of  n  and  is  identical  to  the  solution 
obtained  in  our  previous  one-dimensional  calculations.  We  also  find  from 
Eqs.  (2.1)  and  (2.3) 
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and 


a 


y 


(3.12) 
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a  i 
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Therefore,  for  the  special  case  in  which  the  plasma  conductivity  can  be 
approximated  as  indicated  in  Eq.  (3.10),  the  following  situation  exists.  The 
magnetic  induction  varies  only  in  the  direction  of  propagation  and  no  trans¬ 
verse  currents  exist.  Furthermore,  although  the  conductivity  is  a  function  of 
n,  the  electric-field  intensity  also  depends  on  n  and  in  such  a  way  that  the 
current  density  is  independent  of  n.  Considerable  simplification  results  in 
the  numerical  solution  of  the  equations  if  the  approximation  above  is  made 
and  for  the  problem  at  hand,  as  we  shall  see,  the  approximation  is  not  un¬ 
reasonable.  We  now  turn  to  a  solution  of  the  equations  in  several  cases  of 
interest. 


B.  Limiting-Case  Analytic  Solution 

Although  we  will  be  concerned  primarily  with  the  numerical  solution  of 
the  governing  equations  for  reasonably  general  cases,  it  is  nonetheless  of  in¬ 
terest  to  carry  out  a  limiting-case  analytic  solution.  Such  a  solution  is 
useful  for  checking  the  computational  results  for  the  more  general  cases,  for 
providing  some  physical  insight  into  the  nature  of  the  solution,  and  for  ob¬ 
taining  some  approximate  scaling  relations  which  will  be  discussed  presently. 
In  order  to  make  the  analysis  possible  we  neglect  the  position  dependence  of 
the  plasma  conductivity,  radiation  mean  free  path,  and  degree  of  ionization 
(i.e.,  Xj  and  X2).  Despite  the  apparent  lack  of  justification  for  these  as¬ 
sumptions,  the  analytic  results  will  be  seen  to  be  in  agreement  with  the  more 
general  numerical  results  to  within  at  least  a  few  percent. 

Under  the  assumptions  discussed  above,  we  have  from  Eq.  (3.11) 


B  *  yj  (1-0, 


(3.14) 


and  from  Eq.  (3.9), 


20 


(3.15) 


P  -  wj2  (5-52/2)  . 


The  current  density,  J,  then  follows  directly  from  Eq.  (3.12),  namely. 


j  -  j/ta  . 


(3.16) 


Using  Eqs.  (2.20)  and  (3.16)  in  the  energy  equation,  Eq.  (3.4),  we  find 


2  2  2 

1  +  i  ift*  _j o 
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(3.17) 


where 


t  =  T 


(3.18) 


The  solution  to  Eq.  (3.17)  can  be  written 


/  dz'  /  dn'  G(e,  n;r,n') 


(3.19) 


where  G  is  the  Green  function  for  a  rectangular  surface.  One  hasJ 


G(S,n;e',n')  ■  J 

a  k,p=l 


sin  (kfrg )  sin  (kirg  ^ )  sin  (pun )  sin  (p^n  "  ] 


,.  2  .  2  2.2. 

(k  +  p  /w  ) 


Carrying  out  the  integration,  we  find 


(3.20) 


.2  00  00 


sin  (kir£)  sin  (wrn) 
„  2  .  _2„2  .2. 


T  _  T  j  +  _ £  £  siniKTT^jsin (.pimj 

b  ir4  %aX  T4  k=l , 3, 5  p=i,3,S  kp(k2  +  pV/w2) 

The  mass  density  then  follows  directly  from  Eq.  (2.16)  and  (3.15), 

mQ  M32  (Z  -  C2/2) 
p  =  (1  *  Xj  +  2x2)  kRT  ’ 


(3.21) 


(3.22) 


with  T  given  by  Eq.  (3.21).  In  principle,  Eq.  (3.22)  can  now  be  substituted 
into  Eq.  (2.22)  and  the  resulting  integral  evaluated.  That  result  would  then 
be  used  in  conjunction  with  Eqs.  (3.21)  and  (3.22)  to  determine  the  unknowns 


10.  R.  C our  cent  and  D.  Hilbert ,  Methods  of  Mathematiaal  Phyeice  (Intereoienae, 
New  York ,  1953),  Vol.  I,  Chap.  V. 
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T,  p,  and  £  independently.  In  practice  the  resulting  integral  cannot  be  eval- 

cl 

uated  nor  can  £  be  separated  from  the  infinite  series  in  Eq.  (3.21).  There- 

A 

fore,  a  numerical  solution  is  necessary  even  in  solving  the  equations  resulting 
from  this  analytic  treatment.  In  two  special  cases,  however,  an  exact  analytic 
solution  can  be  found  and  we  now  consider  these  cases. 


1.  Case  I:  vr>>jla  (one-dimensional  limit) 

2  2  2  2 
For  w>>£  ,  p  £  /w  is  negligible  with  respect  to  k  in  Eq.  (3.21) 

for  any  contributing  values  of  p  and  k  and  can  be  neglected.  Furthermore,  one 
can  show  that 


00 
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k=l, 3, 5 


sin 
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(PTtC) 
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*3  (i-s) 
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and 


(3.23) 


j  ii2_iE2l  „  „/4  ,  (3.24) 

P-1,3,5  P 

so  Eq.  (3.21)  becomes 


T  =  j  [i  +  -ij — LLk.Q.1 

b  b  8a s  Xa  T  4j 


1/4 


and  Eq.  (3.22)  becomes 


P  = 


Pj2  (Z-  Z2/2) 


(1+x1+2x2)  kfi  Tb 


1  + 


3j2  5(1-5) 


77 4 


8os  oX  TbJ 


(3.25) 


(3.26) 


The  arc  length  £  can  be  calculated  now  by  substituting  Eq.  (3.26)  into  Eq. 

(2.22)  and  evaluating  the  resulting  integral.  For  purposes  of  carrying  out  the 
integration,  the  first  term  in  brackets  in  Eq.  (3.26)  is  neglected  with  respect 
to  the  second  since  the  latter  term  dominates  except  very  close  to  5  ■  0  and 
5-1.  We  find 


£ 

a 


20  (l+x1+2x2)  kB  Pjl  T  (1/2) 
_ a _ 

13  m0  uj3'2(8os  aX/3)1/4  r2(3/4) 


(3.27) 


2 

where  T  is  the  gamma  function;  the  ratio  r(l/2)/r  (3/4)  is  about  1.2.  There¬ 
fore,  Eqs.  (3.14),  (3.15),  (3.25),  (3.26),  and  (3.27)  provide  uncoupled  expres¬ 
sions  for  B,  p,  T,  p  and  £  . 
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2.  Case  II:  w<<fc 

a 


The  assumption  thht  w<ot  i3  the  one  most  frequently  encountered  in 

8  2 

experiments  of  interest.  Under  this  condition  k  can  be  neglected  with  res¬ 
pect  to  p2S,  2/w“  in  Eq.  (3.21)  and,  using  Eqs.  (3.23)  and  (3.24)  as  before, 
a 

we  find 


T»Tb 
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+  3j2  n  (l-n)w2 

8cc  a  XT.  I 
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(3.28) 


and 


P  * 


n>0  Mj2  (C-52/2) 


(1  +  Xj  +  2x2)  kRTb 


L1  *g  rtT„v  J 
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Similarly,  calculating  the  arc  length  as  before,  we  have 


a  = 
a 


3  (1  +  Xj  +  2x2)  kg  r(i/2)  w1/2 
2m0  u  j3/2  (8as  aX/3)1/4  T2(3/4)  H&1/2 


whence 


(3,29) 


l 


a 


3  (1  +  x1  +  2x2)kB  pz  r Cl/2) 
_ a _ 

2m0  u  r2  (3/4) 


2/3 

1/3 

w _ 

j(8asaX/3)1/6 


(3.30) 


Equation  (3.50)  can  now  be  substituted  into  Eq,  (3.28)  to  produce  T  and  into 
Eq.  (3.29)  to  produce  p. 


3.  Comparison  of  Two  Solutions 

It  is  interesting  to  compare  the  solution  obtained  in  Sec.  B.l  (one¬ 
dimensional  case)  with  that  obtained  in  Sec.  B.2  (two-dimensional  case)  in  the 
light  of  our  remarks  in  Sec.  I  concerning  the  expected  differences  in  the  two 
calculations.  Thus,  we  see  that  in  this  simplified  case  identical  expressions 
are  obtained  for  the  pressure  P  and  magnetic  induction  B.  For  a  constant  arc 
length  l  ,  however,  the  temperature  in  the  interior  of  the  arc  is  smaller  in 

the  two-dimensional  case  than  that  in  the  one-dimensional  case  by  a  factor 
which  decreases  with  increasing  arc  length.  The  latter  result  can  be  seen  by 
comparing  Eqs.  (3.25)  and  (3.28)  in  the  preceding  analysis. 
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C.  Numerical  Technique 

In  the  remainder  of  this  section  we  will  discuss  more  general  solutions 
to  the  governing  equations  for  several  cases.  In  particular,  it  will  no  longer 
be  assumed  that  the  conductivity,  mean  free  path,  and  degree  of  ionization  are 
uniform  in  the  arc  but,  rather,  their  spatial  dependence  will  be  accounted  for, 
Because  of  the  additional  complication  introduced,  we  must  resort  to  numerical 
solutions  of  the  equations,  and  the  particular  numerical  procedure  is  described 
in  this  section. 

We  first  noted  that  the  previous  one- dimensional  results  suggested  that 
substantial  variations  in  the  flow  parameters  could  be  expected  very  near  the 
projectile  surface,  i.e.,  in  the  vicinity  of  5  =  1.  This  problem  was  dealt 
with  in  the  one- dimensional  treatment  by  simply  making  the  step  size  small. 

In  two  dimensions,  however,  such  a  remedy  is  not  practical.  Instead,  we  per¬ 
formed  a  "stretching  transformation"  by  letting 

5  =  1-r2  ,  (3.31) 

and  solved  the  equations  as  a  function  of  5'.  As  can  be  easily  seen, for  con¬ 
stant  &Z'  such  a  transformation  has  the  effect  cf  decreasing  the  step  size  near 
5=1  where  rapid  variations  occur,  and  of  increasing  the  step  size  near  5=0 
where  variations  are  slower. 

We  then  employed  an  iterative  technique  similar  to  that  used  in  our  pre¬ 
vious  one-dimensional  calculations.  Specifically,  we  divided  the  5'n  plane  in¬ 
to  a  grid  and  approximated  all  derivatives  by  finite  differences.  The  result¬ 
ing  coupled  equations  were  then  written  in  the  form 

VR  (i, j )  =  F(Vr  V2  ...Vk'...VN;5'n)  (3.32) 

t*  Vi  *.1% 

where  V^Ci.j)  denotes  the  k  unknown  (P,T,  etc.)  at  the  (i.j)^  point  and 

where  F  denotes  some  function  of  the  remaining  unknowns,  5',  and  n.  The  prime 
on  in  the  argument  of  F  denotes  that  \'k(i,j)  does  not  appear  on  the  right- 

hand  side  of  the  equation.  Initial  values  for  the  variables  on  the  right-hand 
side  of  the  equations  were  obtained  from  the  analytic  approximation  discussed 
in  Sec.  IIIB,  and  new  values  were  calculated  on  the  left-hand  sides.  For  the 
second  iteration,  a  weighted  average  of  the  old  and  new  values  was  used  on  the 
right-hand  sides.  The  process  was  then  repeated  and  was  found  to  converge  in 
a  few  hundred  iterations. 

We  have  performed  calculations  in  which  we  used  the  general  expression  for 

the  plasma  conductivity  in  Eq.  (2.23),  as  well  as  calculations  in  which  we  made 

the  separable-conductivity  approximation  in  Eq.  (3.10).  In  the  latter  case, 

a  and  a  were  obtained  from  the  relations 
x  y 

a  ($)  ~  a  (?»n)  dn  /<o(5,r0>1//2 
x  0 


(3.33) 


(3.34) 


a  (n)  =  /  a(S,n)d£/<otS,Ti>1/’2 
Y  0 

where  o(5,n)  is  given  by  Eq  (2.23)  and  where  the  brackets  denote  an  average 
taken  over  the  entire  arc.  In  all  these  calculations  we  monitored  both  the 
actual  and  the  factored  conductivity  and  good  agreement  was  found  between  the 
two  expressions  in  all  cases.  Presumably,  the  l'eason  for  the  excellent  agree¬ 
ment  is  that  the  conductivity  is  a  very  slowly  varying  function  in  the  inter¬ 
ior  of  the  arc,  and  varies  rapidly  only  with  k  or  n,  but  not  both,  near  the 
boundaries.  The  advantage  of  using  the  factored  conductivity  is  that  the  iter¬ 
ative  procedure  converges  far  more  rapidly  than  in  the  more  general  case,  no 
doubt  because  of  the  absence  of  small  transverse  currents.  In  most  of  the  re¬ 
mainder  of  the  discussion  we  will  confine  ourselves  to  results  in  which  the 
separable-conductivity  approximation  was  made.  In  Sec.  IV,  however,  some  qual¬ 
itative  results  of  the  more  general  case  are  discussed. 

We  have  performed  calculations  using  a  variety  of  step  sizes  for  AC'  and 
An  and  have  found  that  the  values  AC'  -  An  =  0.02  gave  acceptable  results  in 
a  reasonable  calculation  time.  When  AC'  or  An  was  increased  by  a  factor'  of 
two,  for  example,  we  observed  changes  in  some  variables  by  as  much  as  3  or  4% 
from  results  obtained  at  the  smaller  step  size.  Most  quantities  were  in  sub¬ 
stantially  better  agreement  than  that,  however,  and  we  assumed  that  values  of 
AC'  and  An  given  by  0.02  were  sufficiently  small. 

As  has  been  emphasized  before,  the  fluid-dynamic  variables  within  the  arc 
are  nearly  independent  of  provided  the  boundary  temperature  is  small  with 

respect  to  the  maximum  temperature  within  the  arc.  Under  these  conditions 
the  arc  temperature  is  determined  by  the  rate  of  energy  dissipation  only,  and 
the  boundary  temperature  has  negligible  effect.  On  the  other  hand,  in  the  num¬ 
erical  calculations  it  is  desirable  to  have  the  boundary  temperature  as  large 
as  possible  in  order  to  avoid  numerical  difficulties  that  may  arise  from  steep 
temperature  gradients  near  the  boundaries.  Therefore,  in  the  numerical  calcu¬ 
lations  we  set  T^  equal  to  some  fraction  of  the  estimated  maximum  temperature 

within  the  arc,  and  then  lowered  the  value  until  further  reduction  produced  no 
appreciable  change  in  the  arc  properties.  In  practice,  reasonable  values  of 
the  boundary  temperature  were  found  to  be  about  50%  of  the  maximum  temperature 
in  the  arc. 


Analysis  of  Rashleigh-Marshall  Experiment 


In  a  recent  experiment  1  Rashleigh  and  Marshall  (RM)  have  employed  an  arc- 
driven  rail  gun  to  accelerate  a  3-g  projectile  to  about  6  km/s.  The  arc  con¬ 
sisted  of  partially  ionized  copper  vapor.  In  an  effort  to  analyze  the  arc  in 
that  experiment,  we  have  solved  the  two-dimensional  equations  for  determining 
the  arc  properties  for  a  set  of  input  data  roughly  equivalent  to  the  RM  data. 


The  appropriate  data  are  shown  in  Table  I.  All  quantities  were  either 
taken  directly  from  the  experimental  specifications  or  easily  calculated  there¬ 
from.  The  parameter  h  is  an  effective  rail  height  necessary  for  use  in  the 

TT.  S.  C.  Rashleigh  and  R.A.  Marshall,  "Electromagnetic  Acceleration  of  Macro- 
particles  to  High  Velocities J.  Appl.  Phys.  4£,  2540  (1978). 
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TABLE  l.  Experimental  Data  for  Rashle igh-Marshall  Experiment. 


Quantity 

Description 

Value 

w 

Rail  separation 

1.27  x  10‘2  m 

V 

Plasma  height  on  rails 

1.27  x  10~2  m 

hr 

Rail  height 

1.91  x  1Q~2  m 

h 

Effective  rail  height 

1.56  x  10'2  m 

m 

P 

Projectile  mass 

3  x  10‘3  kg 

i 

Pulsed  current 

3  x  105  A 

j 

Current  per  unit  height 
on  rails 

1.92  x  107  A/m 

V0 

Muzzle  voltage 

160  V 

m 

a 

Arc  mass 

2  x  10'4  kg 

mo 

Ion  (or  atom)  mass 

1.1  x  10’25>kg 

a 

Average  acceleration 

6  x  106  m/s2 

P£ 

(See  Eq.  (3.5)) 

0.51  kg/m2 

\ 

CSee  Eq.  (3.5)) 

15.1  kg/m2 

calculations  since  the  quantities  hp  and  hr>  defined  in  the  table,  were  assumed 

to  be  equal  in  the  model.  As  can  be  seen,  h  was  chosen  to  be  the  geometric 

mean  of  the  two  experimental  values.  The  arc  mass  m  was  unknown  experiment- 

Cl 

ally  but  we  used  a  value,  found  by  trial  and  error,  which  gave  an  arc  length 
of  about  10  cm.  Finally,  the  experimental  quantity  Vq  represents  the  poten¬ 
tial  measured  across  the  end  of  the  rails.  The  theoretical  value  can  be  ob¬ 
tained  from  Eq.  (3.13)  by  evaluating  the  expression  at  £  ■  1  and  noting  that 


w  /  E  (l,n)  dn 
0 


(3.35) 


Typical  results  obtained  from  the  ntmerical  calculation  are  shown  in  Table 
II  for  quantities  which  are  independent  of  position.  The  acceleration  a  is 
about  a  factor  of  two  and  one-half  higher  than  the  experimental  value  of  6  x 
ft  2 

10  m/s  .  The  lack  of  agreement  no  doubt  results  because  the  use  of  infinitely 
high  rails  overestimates  the  acceleration  imparted  to  the  projectile,  and 
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TABLE  II.  Results  of  Numerical  Calculations  for  Rashleigh-Marshall 
Experiment 


Quantity 

Value 

a 

1.5  x  107  m/s2 

i 

a 

9.  8  cm 

<Z> 

1.61 

<T> 

3.7  x  104  °K 

A 

V 

c 

V 

1 . 4  x  10  m 

vo 

50  volts 

because  of  the  somewhat  arbitrary  manner  in  which  the  effective  height  h  was 
chosen.  In  a  similar  manner  the  calculated  potential  VQ  is  about  a  •Factor  of 

three  lower  than  the  experimental  value  shown  in  Table  I. 

The  remaining  quantities  in  the  table  are  largely  self  evident.  It  is  of 

interest,  however,  to  compare  them  with  the  corresponding  results  obtained  in 

the  previous  one-dimensional  treatment  for  an  arc  of  roughly  the  same  length. 

As  expected,  we  find  that  the  mean  temperature  of  37, 000° K  is  substantially 

lower  than  for  the  one-dimensional  case  in  which  we  found  <T>  -  57,000°.  As  a 

result  of  the  lower  temperature,  the  mean  ionization  <Z>  is  also  lower  than 

the  one-dimensional  value  of  1.93.  Since  the  pressure  should  be  nearly  the 

same  in  the  two  calculations,  as  we  indicated  before,  the  mass  density  in  the 

arc  should  be  higher  than  in  the  one-dimensional  case  in  order  to  satisfy  the 

equation  of  state  and,  therefore,  the  electron  density  is  also  higher.  Thus, 

25-3  .  . 

we  found  <n  >  =  9.9  x  10  m  in  our  previous  treatment  and  this  value  is  m- 

®  26  -  3 

deed  smaller  than  the  present  value  of  1.4  x  10  m”  .  Finally,  we  observe 

that  the  potential  is  roughly  the  same  in  the  two  calculations  (50  volts 

compared  to  47  volts).  It  is  perhaps  surprising  that  the  potential  is  not  sig¬ 
nificantly  higher  than  that  in  the  one-dimensional  case  since  the  conductivity 

3/2 

of  the  arc  varies  as  T  .  However,  in  this  particular  temperature  range,  the 
degree  of  ionization  is  also  smaller  and  there  are,  consequently,  fewer  colli¬ 
sions  among  the  electrons.  The  reduced  scattering  tends  to  increase  the  con¬ 
ductivity  and  this  effect  nearly  compensates  for  the  effect  of  the  reduced 
temperature.  Thus,  little  change  in  the  conductivity  results. 

Graphs  of  the  position-dependent  quantities  are  shown  in  Figs.  2-6.  For 
the  pressure  and  magnetic  induction,  graphs  are  drawn  as  a  function  of  £. 

These  quantities  do  not  depend  on  n  when  the  conductivity  is  approximated  by 
a  separable  function,  and  only  weakly  so  when  the  more  general  expression  is 
used.  This  behavior,  incidentally,  is  indicative  of  our  earlier  statement 
that  only  small  differences  between  these  parameters  should  be  expected  in  the 
one-  and  two-dimensional  calculations. 


The  pressure,  shown  in  Fig.  2,  rises  monotonically  from  zero  at  the  back 
of  the  arc  to  about  220  MPa  at  the  projectile  surface.  Its  functional  depend¬ 
ence  may  be  compared  with  the  approximate  analytic  form  in  Eq.  (3. IS).  Simi¬ 
larly,  the  magnetic  induction  field  varies  from  about  24  Tesla  at  the  trailing 
edge  of  the  arc  to  zero  at  the  projectile  surface.  If  the  conductivity  of  the 
arc  were  constant,  the  variation  between  these  two  values  would  be  linear. 

The  fairly  weak  dependence  of  the  conductivity  on  position,  however,  gives  rise 
to  the  small  deviations  from  linear  behavior  observed  in  the  graph. 

In  Fig.  4,  two-dimensional  effects  of  the  temperature  profile  are  shown 
by  plotting  isotherms  within  the  arc.  Here,  unlike  for  the  pressure  and  fields, 
significant  variations  of  the  temperature  were  observed  as  a  function  of  n,  par¬ 
ticularly  near  the  surface  of  the  rails.  It  may  furthermore  be  noted  that  there 
is  a  high  concentration  of  energy  near  the  projectile  surface,  with  a  very  steep 
gradient  as  the  surface  is  approached.  The  high  temperature  near  the  surface 
can  be  explained  by  noting  that  there  is  an  increase  in  mass  density  with  in¬ 
creasing  K  (see  Fig.  5)  and,  consequently,  a  decrease  in  the  radiation  mean 
free  path.  Thus,  photons  created  in  the  right-most  part  of  the  arc  are  radiated 
away  with  more  difficulty  than  those  in  the  left -most  part.  Once  the  photons 
are  within  a  mean  free  path  or  so  of  the  boundary,  however,  they  can  essenti¬ 
ally  "see''  the  boundary  and  are  easily  radiated  outward.  Therefore,  the  tem¬ 
perature  drops  rapidly  very  near  the  projectile  surface.  The  maximum  tempera¬ 
ture  in  the  arc  was  found  to  be  about  46,0008K  whereas  in  our  previous  one¬ 
dimensional  work  the  maximum  was  observed  to  be  about  70,000  deg.  As  noted 
earlier,  the  temperature  was  expected  to  be  smaller  for  the  two-dimensional 
calculation. 

In  Fig.  5  are  plotted  lines  of  constant  mass  density.  The  density  may  be 
seen  to  vanish  at  the  back  of  the  arc  and,  as  in  the  temperature  plot,  a  rather 
rapid  variation  in  the  density  is  observed  near  the  projectile  surface.  The  in¬ 
crease  there  results  because,  as  the  temperature  falls,  the  degree  of  ionization 
also  falls,  and  both  these  effects  tend  to  reduce  the  pressure  at  the  projectile 
surface.  Therefore,  a  steady  pressure  can  be  maintained  at  the  surface  only  if 
there  is  a  proportionate  increase  in  density  there. 

Finally,  lines  of  constant  electron  density  may  be  observed  in  Fig.  6. 

The  curves  are  similar  to  those  for  the  mass  density  but  the  gradient  near 
5  =  1  is  significantly  less  pronounced.  The  reason  for  the  smaller  gradient, 
of  course,  is  that  the  ions  near  the  projectile  surface  are  less  highly  ionized 
than  those  in  the  interior  of  the  arc  because  of  the  lower  temperature. 

E.  Analysis  of  Westinghouse  Experiment 

In  addition  to  the  analysis  of  the  arc  in  the  Rashleigh-MarshaJ 1  experi¬ 
ment,  we  have  also  performed  calculations  to  determine  the  properties  of  an 
arc  for  a  rail  gun  comparable  in  size  to  that  under  construction  at  Westing- 
house.  Typical  input  parameters,  again  taken  largely  from  the  experimental 
specifications,  are  shewn  in  Table  III.  A?  in  the  RM  case,  the  parameter 

was  chosen  by  trial  and  error  to  produce  an  arc  length  of  roughly  ten  centi¬ 
meters;  the  pulsed  current  i  was  chosen  so  as  to  produce  an  acceleration  suf¬ 
ficient  for  the  projectile  to  reach  a  velocity  of  3  km/s  in  a  distance  of  4m. 


TABLE  III.  Experimental  Data  for  Westinghouse  Gun. 


Quantity  Description  Value 


K 

Rail  separation 

5.0  x  10‘2 

m 

h 

Rail  height 

5.0  x  10"2 

m 

m 

P 

Projectile  mass 

0.3  kg 

i 

Pulsed  current 

7.35  x  105 

A 

j 

Current  per  unit  height  on  rails 

1.47  x  107 

A/m 

m 

a 

Arc  mass 

1.2  x  10“3 

kg 

P£ 

(See  Eq.  (3.5)) 

0.48  kg/m2 

a 

(See  Eq.  (3,5)) 

120  kg/m2 

Position-dependent  quantities  in  the  calculation  vary  in  essentially  the 
same  manner  as  for  the  RM  case,  so  the  significant  results  will  be  presented 
only  in  tabular  form.  These  results  are  shown  in  Table  IV.  The  results 

TABLE  IV.  Results  of  Numerical  Solution  for  Westinghouse  Gun. 


Quantity  Value 


a 

1.13  x 

:  106  m/s 

n 

a 

9.9  cm 

o 

> 

171  volts 

<T> 

4.2  x 

O 

O 

T 

max 

5.2  x 

US 

o 

O 

H 

<n  > 
e 

8.0  x 

in25  "3 
10  m 

ne  max 

1.6  x 

. n26  -3 

10  m 

p 

max 

17.0  kg/m 

<Z> 

1.83 

P 

max 

1.4  x 

108  Pa 

B 

max 

18.4  T 

should  be  compared  with  those  in  both  Table  II  and  in  Figs.  2-6.  The  magni¬ 
tudes  of  the  various  flow  parameters  are  sufficiently  similar  to  those  obtained 
in  the  RM  experiment  to  suggest  that  no  major  difficulties  should  be  expected 
in  using  an  arc  armature  in  this  larger  gun.  It  is  Interesting  to  note,  how¬ 
ever,  that  the  two-dimensional  model  predicts  a  slightly  higher  temperature  for 
the  large-gun  arc  than  for  the  RM  arc,  the  opposite  of  the  prediction  in  the 
one- dimensional  case.  The  reasons  for  the  contradictory  predictions  as  well 
as  some  other  results  of  the  calculation  will  be  made  clearer  in  the  following 
section  where  some  approximate  scaling  laws  will  be  derived. 

F.  Scaling  Relations  and  Optimum  Arc  Materials 

The  numerical  calculations  undertaken  previously  are  difficult  to  carry 
out  and  computationally  rather  inefficient.  It  is  of  interest,  therefore,  to 
attempt  to  derive  some  approximate  scaling  relations  which  indicate  how  arc 
characteristics  vary  with  experimental  parameters.  Such  a  set  of  relations  can 
be  derived  from  the  special-case  analytic  solution  presented  in  Sec.  IIIB  if 
some  additional  assumptions  are  made,  and  if  either  Z  <<w  or  z  >>w.  In  par- 

ticular,  we  make  the  following  additional  approximations:  (1)  In  order  to  ac¬ 
count  to  some  extent  for  the  temperature  dependence  of  the  conductivity  in  the 
results  of  Secs.  IIIB.l  and  IIIB. 2,  we  assume 

a  a  T‘I2I1  (3.36) 


and  neglect  the  logarithmic  term  in  Eq.  (2.23).  (2)  We  note  that  the  radiation 

mean  free  path  X  in  Eq.  (2.21)  varies  approximately  according  to  the  relation6 


X  a 


2  T7/2 

m0  T 


2  . 

P  <J> 


(<j>  +  1)' 


(3.37) 


where  <j>  for  this  case  can  be  approximated  by 


<j>  =  Xj  +  2x2 


(3.38) 


(3)  Finally,  we  note  that  the  parameter  <j>  can  be  well  approximated  by  the 
parameter  Z  (see  Eq.  (2.24)  if  either  single  or  double  ionization  is  dominant. 
In  order  to  simplify  the  results  as  much  as  possible,  we  will  assume  this  to 
be  the  case. 

If  we  make  use  of  the  assumptions  noted  above  in  the  results  of  Secs. 
IIIB.l  and  IIIB. 2  and  solve  the  coupled  equations  for  the  unknown  quantities, 
we  obtain  the  scaling  relations  indicated  in  Table  V.  The  results  should  be 
reasonably  valid  provided  they  are  applied  not  too  close  to  the  arc's  boundar¬ 
ies.  In  order  to  use  the  table  to  predict  some  quantity  of  interest  for  a  case 
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TABLE  V.  General  Scaling  Factors  . 
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in  which  no  numerical  calculations  exist,  the  following  procedure  is  followed. 
The  value  of  a  particular  quantity  in  the  RM  experiment  (column  1)  is  multi¬ 
plied  by  ratio  of  the  scaling  factors  for  the  unknown  and  RM  cases.  The  re¬ 
sult,  then,  is  the  value  of  the  parameter  of  interest  in  the  unknown  case. 

The  degree  of  ionization,  Z,  is  not  easily  expressible  as  a  function  of  the 
experimental  quantities  as  are  the  other  unknowns.  Fortunately,  however, 
large  changes  in  this  parameter  do  not  occur  even  with  fairly  significant 
changes  in  the  experimental  conditions,  and  it  can  often  be  treated  as  a  con¬ 
stant  in  the  two  cases  of  interest. 


The  results  in  Table  V  are  general  and,  except  for  their  dependence  on  Z, 
depend  only  on  the  experimental  conditions,  In  the  calculations  undertaken  so 
far,  however,  we  have  varied  the  experimental  parameter,  ,  in  order  to  main- 

a 

tain  a  constant  arc  length  of  roughly  10  cm.  It  is,  therefore,  of  interest  to 

determine  scaling  relations  as  a  function  of  arc  length  and  treat  l  as  an  ex- 

& 

perimentally  determined  quantity.  To  do  so,  we  solve  for  p  in  terms  of  £ 

a 

(see  the  sixth  entry  in  Table  V)  and  substitute  the  result  elsewhere  in  the 
table.  We  then  obtain  the  results  shown  in  Table  VI. 


TABLE  VI.  Scaling  factors  as  a  Function  of  Arc  Length. 
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To  demonstrate  the  application  of  the  table  we  have  used  it  to  determine 
approximate  results  for  the  Westinghouse  experiment,  using  the  numerical  re¬ 
sults  for  the  RM  case.  Shown  in  column  2  of  Table  VII  are  numerical  results 
for  the  RM  experiment.  In  column  three  we  have  repeated,  for  convenience,  the 
numerical  results  for  the  Westinghouse  case  (see  Table  IV).  In  column  four 
are  numerical  values  of  the  scaling-factor  ratios  (see  Tables  I,  III,  and  vi) . 
Finally,  in  column  five  are  the  results  for  the  Westinghouse  experiment  as  pre¬ 
dicted  by  the  scaling  relations.  Comparing  columns  three  and  five  we  see  that, 
except  for  the  potential  Vn  and  the  electron  density  n  ,  the  agreement  is  ex- 

cellent.  The  lack  of  agreement  in  the  last  two  quantities  results  because 
they  are  fairly  sensitive  functions  of  Z,  and  this  parameter  does  change  some¬ 
what  in  the  two  cases  (see  Tables  II  and  IV  for  comparison  of  actual  numerical 
values).  In  addition,  the  condition  2.  >>w,  required  for  the  validity  of  these 

relations,  is  not  strictly  valid  for  the  case  at  hand. 
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Finally,  it  is  of  interest  to  ask  what  the  scaling  relations  predict  con¬ 
cerning  the  optimum  materials  to  use  for  the  arc.  In  particular,  it  has  been 
suggested  in  the  past  that  it  might  be  desirable  to  use  some  substance  that 
is  more  easily  ionized  than  copper,  say  potassium  or  cesium,  and  attempt  to 
obtain  equivalent,  conductivity  at  substantially  lower  temperatures.  Unfor¬ 
tunately,  as  can  be  seen  from  Table  VI,  use  of  these  substances  having  low 
ionization  potentials  will  probably  not  make  any  appreciable  difference  in  the 
final  temperature  of  the  arc.  The  reason  is  that  temperatures  associated  with 
the  arc  in  any  practical  experiment  are  always  sufficiently  high  that  almost 
complete  ionization  of  the  arc  is  obtained.  Once  the  condition  has  been 
achieved,  additional  increases  in  the  electron  density,  which  would  accompany 
lower  ionization  potential,  do  not  result  in  higher  conductivity  and  lower  tem¬ 
peratures.  In  fact,  as  can  be  seen  from  I'q.  (2.23),  as  2  becomes  larger,  the 
conductivity  actually  decreases.  Referring  to  Table  VI,  it  is  evident  that 
the  temperature  is  controlled  largely  by  the  current  and,  to  a  significantly 
smaller  extent,  by  the  length  of  the  arc,  the  rail  separation,  and  the  degree 
of  ionization  Z.  High  currents  are,  of  course,  necessary  to  achieve  the  high 
accelerations  required  so  we  have  little  control  over  this  parameter.  Very 
small  changes  in  the  temperature  cart  be  obtained  by  varying  the  other  three 
parameters,  however,  and  it  is  generally  desirable  to  make  the  length  of  the 
arc  as  long  as  possible,  the  rail  separation  as  small  as  possible,  and  to 
choose  as  the  arc  material  a  substance  having  a  high  second- ionization  poten¬ 
tial.  Again,  however,  making  the  optimum  choices  is  not  going  to  change  the 
temperature  significantly  as  is  evident  from  the  table.  Doubling  the  arc 
length,  for  instance,  for  a  situation  in  which  Z  is  held  constant, results  in  a 
temperature  reduction  of  only  about  12°6. 

To  demonstrate  the  validity  of  these  remarks,  v/e  have  also  performed  some 
numerical  calculations  for  the  RM  gun  using  potassium  vapor  for  the  arc,  rather 
than  the  copper  vapor  used  previously.  The  first  and  second  ionization  poten¬ 
tials  for  potassium  are  4.3  eV  and  31.8  eV,  respectively,  and  its  atomic  mass 
is  39.  Corresponding  values  for  copper  are  7.7  eV,  20.3  eV,  and  64.  For  con¬ 
stant  arc  length  10  cm),  very  small  differences  were  noted  in  the  tempera¬ 
ture  in  the  two  calculations.  In  fact,  the  only  parameter  in  Table  VI  which 
changed  more  than  a  few  percent  was  the  mass  density  p.  It  was  found  to  be 
about  a  factor  of  two  smaller  for  the  potassium  arc  than  for  the  copper  arc. 

The  result  is  predictable  from  Table  VI  and  can  be  ascribed  to  the  smaller 
atomic  mass,  m^,  for  potassium. 


IV.  DISCUSSION 

In  Ref.  1  we  discussed  in  some  detail  the  assumptions  employed  in  the  one¬ 
dimensional  calculation.  For  the  most  part,  that  discussion  applies  to  the 
two-dimensional  case  as  well.  Therefore,  in  this  section  we  will  confine  our¬ 
selves  to  summarizing  the  results  of  the  two-dimensional  calculation  and  to 
discussing  only  assumptions  pertinent  to  that  case. 

We  have  derived  a  general,  time-dependent  formalism  for  studying  arc  dyna¬ 
mics  in  two  dimensions.  For  purposes  of  carrying  out  the  calculation,  we  have 
neglected  the  effect  of  the  arc's  acceleration  upon  its  fluid-dynamical  charac¬ 
teristics  »  and  have  assumed  that  those  characteristics  are  steady  in  a  frame 
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of  reference  moving  with  the  arc.  The  latter  assumption  was  also  made  in  our 
previous  one-dimensional  treatment.  The  resulting  governing  equations  were 
then  solved  to  analyze  the  arcs  both  in  the  RM  experiment  and  for  a  gun  com¬ 
parable  in  size  to  that  under  construction  at  Westinghouse.  From  approximate 
analytic  solutions  to  the  governing  equations,  some  approximate  scaling  rela¬ 
tions  were  derived  and  these  relations  were  shown  to  yield  results  in  reason¬ 
ably  good  agreement  with  the  more  general  numerical  solutions. 

The  assumption  that  a  is  negligible  in  the  governing  equations  may  be 
viewed  as  the  lowest-order  approximation  in  any  of  an  infinite  number  of 
perturbation-expansion  solutions,  valid  for  small  values  of  the  acceleration. 
For  example,  if  one  assumes  that  all  the  time-derivative  terms  in  the  equations 
(including  the  velocities)  are  small,  as  in  a  quasi-static  approximation,  a 
reasonably  straightforward  expansion  of  the  equations  can  be  obtained,  with  the 
equations  actually  solved  corresponding  to  lowest  order.  Higher-order  approxi¬ 
mations  can  also  be  written  down  in  such  an  expansion,  but  the  results  quickly 

12 

become  quite  tedious.  At  any  rate,  it  is  clear  that  there  is  no  steady  solu¬ 
tion  in  higher  orders  of  the  two-dimensional  problem  for  a  situation  in  which 
t.he  gas  is  at  rest  relative  to  the  projectile.  This  situation  may  be  contrast¬ 
ed  with  the  one-dimensional  case  where  such  an  expansion  may  also  be  carried 
out.  In  this  instance,  however,  one  finds  a  steady  solution  to  the  governing 

12 

equations  in  all  the  higher  orders.  Moreover,  it  is  also  found  in  this  case 
that  the  lowest-order  terms  are  an  excellent  approximation  to  the  exact  solu¬ 
tions;  the  two  results  differ  by  only  a  few  percent. 

We  should  also  point  out  that  a  steady  solution  may  exist  in  the  two- 
dimensional  case  for  a  situation  in  which  the  arc  is  in  motion  relative  to  the 
projectile,  provided  that  motion  is  itself  steady.  One  might  envision,  for  ex¬ 
ample,  two  counter-rotating  flows  in  the  arc  with  the  velocity  as  a  function 
of  position  remaining  constant  in  time.  The  counter-rotation  is  necessary,  of 
course,  to  conserve  the  arc's  total  angular  momentum. 


It  is  of  some  interest  to  attempt  to  estimate  the  conditions  under  which 
the  velocity  might  be  expected  to  be  small  in  such  a  situation.  To  do  so,  we 
first  of  all  note  that  p  ^  p^  /l,  where  L  is  some  distance  characteristic  of 

a 


the  arc  dimensions,  whereas  P  and  the  field-dependent  terms  in  Eqs.  (2.13)  and 

(2.14)  are  of  the  order  of  p^  a.  Consequently,  one  can  neglect  the  term  p^  w a 

a  a 

as  well  as  terms  containing  v'  and  v',  in  the  momentum-conservation  equations 

x  y 

of  interest. 

velocities 
assume 
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<<  1.  This  condition  is  satisfied  in  most  cases 


provided  p  .p 

a  P 

In  order  to  determine  the  conditions  under  which  the  "rotational” 
v  ',  and  v'  are  negligible  in  the  energy  equation,  Eq.  (2.15),  we 
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where  AT  represents  the  change  in  T  in  a  distance  comparable  to  L.  Using 
these  assumptions,  along  with  Eqs.  (2.16)  and  (2.20),  we  estimate  from  Eq. 
(2.15)  that 


(4.2) 


where  v„  stands  for  rotational  velocity.  Clearly,  to  neglect  v we  must  have 
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Finally,  to  determine  conditions  under  which  we  can  neglect  rotational  veloci 
ties  for  purposes  of  determining  the  magnetic-induction  field,  we  assume 
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and  conclude  from  Eq.  (2.6)  that  we  must  have 


(4.5) 


Unfortunately,  it  is  difficult  to  obtain  reliable  estimates  of  v  .  An  upper 
bound  is  probably  some  number  of  the  order  of  r 
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for  this  value  of  v Eqs.  (4.3)  and  (4.5)  are  not  well  satisfied  for  all  cases 

of  experimental  interest.  It  is  likely,  however,  that  v^  is  sizably  smaller 

than  the  value  suggested  in  Eq.  (4.6)  and,  for  more  realistic  values  of  this 
parameter,  Eqs.  (4.3)  and  (4.5)  may  be  easily  satisfied.  Some  further  study  of 
this  problem,  as  well  as  the  question  of  whether  a  steady  state  can  even  exist, 
is  evidently  desirable  and  is  deferred  to  future  work. 
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As  pointed  out  previously,  all  calculations  discussed  here  are  ones  for 
which  the  separable-conductivity  approximation  in  Eq.  (3.10),  was  made.  How¬ 
ever,  we  have  also  carried  out  calculations  using  the  more  general  expression 
for  the  conductivity.  For  these  latter  calculations,  we  obtained  nonzero 
values  of  Jx  as  well  as  some  n  dependence  in  both  B  and  P.  The  magnitude  of 

J^,  however,  was  some  four  orders  of  magnitude  smaller  than  that  of  J  and 

the  x)  dependence  of  both  B  and  P  was  negligible.  These  results  clearly  demon¬ 
strate  that  the  separable-conductivity  approximation  is  quite  good  for  the 
problem  at  hand;  monitoring  both  the  factored  and  the  actual  conductivity  in 
our  numerical  calculations  led  us  to  the  same  conclusion.  As  a  point  of 
academic  interest,  it  is  perhaps  worth  mentioning  that  the  small  values  of  Jx 

are  such  that  they  give  rise  to  a  total  current  density  7  which  flows  toward 
regions  where  the  conductivity  is  high. 

In  Ref.  1  we  pointed  out  that  a  one-dimensional  model  could  be  expected 
to  substantially  overestimate  the  mean  temperature  within  the  arc,  because 
only  two  surfaces  of  the  arc  were  available  for  radiation  in  such  a  model.  It 

was  pointed  out  that  the  mean  temperature  should  vary  roughly  as  A~*^4,  where 
A  is  the  surface  area  available  for  radiation.  Consequently,  we  would  expect 


T3 

sr-  'v 

T1 

A.  1/4 

A3 

(4.7) 

where 

A1  = 

2  h  w 

(4.8) 

A3  " 

2(1  vi  +  l  h  +  hw) 
fl  ft 

(4.9) 

and  where  T^  and  T^  denote  the  mean  temperatures  in  the  one-  and  three- 

dimensional  cases,  respectively.  For  the  RM  data,  we  have  (Aj/A^)1/4  ^  1/2,  so 

one  expects  a  three-dimensional  treatment  to  yield  a  temperature  a  factor  of  two 
or  so  smaller  than  obtained  in  the  one-dimensional  case.  However,  we  also  have 


*  2(hw  +  £a  h)  (4.10) 

1/4 

and  (Aj/Aj,)  v  1.15  for  the  RM  case.  Therefore,  the  expected  difference  in 

temperature  obtained  in  a  two-  and  three-dimensional  treatment  is  only  about 
15%.  Thus,  the  results  of  the  two-dimensional  calculation  are  probably  a  far 
more  accurate  representation  of  the  arc's  actual  properties  than  are  the  one¬ 
dimensional  results  and  they  are,  furthermore,  probably  not  too  far  from  re 
suits  that  would  be  obtained  in  a  three-dimensional  treatment.  It  is  of  in¬ 
terest  to  note  that  for  the  RM  data 
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whereas  the  ratio  of  the  mean  temperature  in  the  two-dimensional  and  one¬ 
dimensional  calculations  (37,000  °K  and  57,000  *K,  respectively)  is  0.65 
The  two  numbers  are  in  fairly  close  agreement. 

In  the  near  future  a  number  of  calculations  will  be  of  interest.  It  is 
first  of  all  desirable  to  consider  higher-order  terms  in  the  governing  equa¬ 
tions  in  order  to  include  the  effects  of  the  rotational  motion  of  the  arc  dis¬ 
cussed  previously.  It  is  not  likely  that  these  effects  will  produce  any  major 
change  in  the  arc's  flu id -dynamic  properties,  but  the  importance  of  the  ef¬ 
fects  does  need  to  be  assessed.  Second,  it  will  be  of  interest  to  study  the 
effects  of  the  plasma  upon  the  projectile  and  rails.  These  calculations  will 
provide  a  first  step  in  evaluating  the  erosive  effects  of  arc.  Finally,  some 
future  consideration  should  be  given  to  time- dependent  effects.  Only  through 
such  treatments  can  we  study  the  formation  and  early  development  of  the  arc 
and  assess  the  validity  of  the  steady-state  assumption.  These  problems  are 
slated  for  study  in  the  near  future. 

NOTE  ADDED  IN  PROOF:  Since  completion  of  this  work,  the  Westinghouse  rail  gun 

13 

has  been  completed  and  fired  several  times.  In  the  most  successful  firing, 
a  317-g  projectile  was  accelerated  to  a  velocity  of  4.2  km/s  in  a  distance  of 
about  5m.  The  resulting  projectile  energy  is  the  largest  yet  achieved  by  a 
rail  gun,  and  the  experiment  represents  a  remarkable  advance  in  rail-gun  tech¬ 
nology. 
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Dr.  Raymond  P.  Gogol ewski 
Dr.  Harry  FaiT 
1400  Wilson  Blvd 
Arlington,  VA  22209 

1  Office  of  Assistant  Secretary  of 
the  Army 

ATTN:  RDA,  Dr.  Joseph  Yang 
Rocm  2E672,  Pentagon 
Washington,  DC  20310 

1  HQDA  (OAMA-ARZ-A,  Dr.  Marvin 
Lasser) 

Washington,  DC  20310 
1  Commander: 

US  Army  Materiel  Development  8 
Readiness  Command 
ATTN:  DRCDMD-ST 
5001  Eisenhower  Avenue 
Alexandria,  VA  22333 

1  Commander 

US  Army  Amrament  Research  and 
Development  Command 
ATTN :  DRDAR-TDC,  Dr.  D.  Gyorog 
Dover,  NJ  07801 


No.  of 

Copies  Organization 
1  Commander 

US  Army  Materiel  Development  <5 
Readiness  Command 
ATTN:  DRCLDC,  Mr.  Lang worthy 
5001  Eisenhower  Avenue 
Alexandria,  VA  22333 

6  Commander 

US  Army  Armament  Research  and 
Development  Command 
ATTN:  DRDAR-TSS  (2  cys) 

DRDAR-LCA,  Mr.  J.A.  Bennett 
DRDAR-LCA,  Dr.  T.  Gora 
DRDAR-LCA,  Dr.  P.  Kenrny 
DRDAR-SCA-W,  R.  Goldstein 
Dover,  NJ  07801 

1  Commander 

US  Army  Armament  Materiel  Readiness 
Command 

ATTN:  DRSAR-LEP-L,  Tech  Library 
Rock  Island,  1L  61299 

1  Director 

US  Army  ARRADCOM 
Benet  Weapons  Laboratory 
ATTN:  DRDAR-ICB-TL 
Watervliet,  NY  12189 

1  Commander 

US  Army  Aviation  Research  § 
Development  Command 
ATTN:  DRDAV-E 
4300  Goodfeilow  Blvd. 

St.  Louis,  M0  63120 

1  Director 

US  Army  Air  Mobility  Research  § 
Development  Laboratory 
Ames  Research  Center 
Moffett  Field,  CA  94035 

1  Commander 

US  Army  Communications  Research  $ 
Development  Command 
ATTN’:  DRDC0-PPA-SA 
Ft.  Monmouth,  NJ  07703 
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1  Commander  4 

US  Army  Electronics  R$D  Command 
Technical  Support  Activity 

ATTN:  DELSD-L 

Tt.  Monmouth,  NJ  07703 

2  Commander 

US  Army  Missile  Command 
ATTN:  DRSMI-R  j 

DRSMI-YDL 

Redstone  Arsenal,  AL  35898 
1  Commander 

US  Army  Tank  Automotive  Research 
6  Development  Command  1 

ATTN:  DRDTA-UI, 

Warren,  MI  48090 


1  Director 

US  Army  TPADOC  Systems  Analysis  2 

Activity 

ATTN:  ATAA-SL,  ^ech  Library 
White  Sands  Missile  Range, 

NM  88002  ] 

2  Commander 

US  Army  Research  Office  1 

ATTN:  Dr.  Fred  Schmiedechoff 
Dr.  M.  Ciftan 

PO  Box  12211  1 

Research  Triangle  Park,  NC  27709 

2  Commander  1 

Naval  Air  Systems  Command 
ATTN:  Dr.  R.J.  Wasneski,  Code  350F 
John  A.  Reif,  AIR  350D  1 

Washington,  DC  20360 


1 


Commander 

Naval  Surface  Weapons  Center 
ATTN:  Henry  B.  Odom,  Code  P-12 

Dr.  M.  Franklin  Rose,  Code  F-04 
P.T.  Adams  5  D.L.  Brunson, 

Code  G-35 

Dahlgren,  VA  22448 
Commander 

US  Naval  Research  Laboratory 
ATTN:  Dick  Ford,  Cede  4774 
Washington,  DC  20301 

HQ  AFSC/XRB 

ATTN:  CPT  Dermis  Kiri in 
Andrews  AF  Base,  MD  20331 


AFATL/DI.DB  CLanny  Burdge) 

(Bill  Lucas) 

Eglin  AF  Base,  FL  32542 

AFATL  (Richard  Walley) 

Eglin  AF  Base,  FL  32542 

AFWL  (Dr.  William  L.  Baker) 

Kirtlund  AFB,  NM  87117 

AFWL/NTYP  (John  Generosa) 

Xirtland  AFB,  NM  87117 

AFAPL  (Dr.  Charles  E.  Oberly) 
Wright- Patterson  AF  Base,  OH  45433 

AFWAL/P00S-2  (CPT  Jerry  Clark) 
Wright  Patterson  AF  Base,  OH  45433 

Director 

Brookhaven  National  Laboratory 
ATTN:  Dr.  James  R.  Powell 
Bldg  129 

Upton,  NY  11973 
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2  Director 

Los  Alamos  Scientific  Laboratory 
ATTN:  Dr.  Clarence  M.  Fowler, 
MS970 

Dr.  Denis  R.  Peterson, 
MS985 

Los  Alamos,  NM  87544 

3  NASA-LeRC  MS  501- 7 
ATTN:  Bill  Kerslake 

Frank  Terdan 
Mike  Brasher 
2100  Brook  Park  RD 
Cleveland,  OH  4413S 

1  International  Applied  Physics 
ATTN:  Dr.  John  P.  Barber 
2400  Glenheath  Drive 
Kettering,  OH  45400 

1  JAYC0R 

ATTN:  Dr.  Derek  Tidman 
205  S.  Whiting  Street 
Alexandria,  VA  22304 

2  R8D  Associates 

ATTN:  Mr.  Ronald  Cunningham 
Dr.  Peter  Turchi 
P0  Box  9695 

Marina  Del  Rey,  CA  90291 


1  Science  Applications,  Inc. 
ATTN:  Dr.  Jad  H.  Batteh 
6600  Powers  Ferry  RD 
Suite  220 

Atlanta,  GA  30339 

1  Science  Applications,  Inc. 
Corporate  Headquarters 
ATTN:  Dr.  Frank  Chilton 
1200  Prospect  Street 
La  Jolla,  CA  92038 


1  Westinghouse  Research  6  Development 

Laboratory 

ATTN:  Dr.  Ian  R.  McNab 
1310  Beulah  RD 
Pittsburgh,  PA  15253 

2  Massachusetts  Institute  of  Technology 
Francis  Bitter  National  Magnet 

Laboratory 

ATTN:  Dr.  Henry  H.  Kolm, 

Mr.  Peter  Mongau 
NW-14-3102 
170  Albany  Street 
Cambridge,  MA  02139 

2  University  of  Florida 

ATTN:  Prof.  K.  T.  Millsaps 
Prof.  B.  M.  Leadon 
Gainesville,  FL  32601 


2  University  of  Texas 

Center  of  Electromechanics 
ATTN:  Dr.  Richard  A.  Marshall 
Mr.  William  F.  Weldon 
167  Taylor  Hall 
Austin,  TX  78712 

Aberdeen  Proving  Ground 

Dir,  USAMSAA 
ATTN:  DRXSY-D 

DRXSY-MP,  H.  Cohen 
Cdr,  USATECOM 
ATTN:  DRSTE-TO-F 
Dir,  US ACS L,  Bldg  E3516,  EA 
ATTN:  DRDAR-CLB-PA 
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